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ABSTRACT 

We use a 64ft,~ 1 Mpc dark matter (DM) only cosmological simulation to examine the large 
scale orientation of haloes and substructures with respect the cosmic web. A web classification 
scheme based on the velocity shear tensor is used to assign to each halo in the simulation a 
web type: knot, filament, sheet or void. Using ~ 10 6 haloes that span ^3 orders of magnitude 
in mass the orientation of the halo's spin and the orbital angular momentum of subhaloes with 
respect to the eigenvectors of the shear tensor is examined. We find that the orbital angular 
momentum of subhaloes tends to align with the intermediate eigenvector of the velocity shear 
tensor for all haloes in knots, filaments and sheets. This result indicates that the kinematics 
of substructures located deep within the virialized regions of a halo is determined by its infall 
which in turn is determined by the large scale velocity shear, a surprising result given the 
virilaized nature of haloes. The non-random nature of subhalo accretion is thus imprinted on 
the angular momentum measured at z = 0. We also find that haloes' spin axis is aligned with 
the third eigenvector of the velocity shear tensor in filaments and sheets: the halo spin axis 
points along filaments and lies in the plane of cosmic sheets. 



1 INTRODUCTION 



Satellite galaxies represent a particularly fascinating subset of the 
low redshift universe. Among their interesting properties is the so- 
called "Holmberg effect" (Holmberg 1969| |Zaritsky et al.||l997J 
which finds that satellites tend to be oriented close to the poles 
of their host disc galaxies, avoiding the planar regions. This type 
of polar alignment is also seen in the Milky Way lKroup a~et al.| 
2005 Metz etal. 2009 1 and was probably first noted over 3 decades 
ago (|Lynden-Bell 19821, an adequate explanation of which is still 
lacking. 

External galaxies (for example in the SDSS) often exhibit a 
spatial alignment of satellites with the long axis of t he host (|Sales| 
' & Lambas||2004{ |Brainerd||20031 |Yang et al.||2006{ |Agustsson &| 
Brainerd 2010) although this alignment is strongest for the older 



redder satellites - younger blue satellite galaxies exhibit no signifi- 
cant alignment. Its not clear whether the Milky Way's polar align- 
ment is at odds with these results, as the surface brightness of many 
of the Milky Way's ultra faint satellite galaxies is many orders of 
magnitude lower than what can be observed extra-galactically. 

In addition to their anisotropic spatial distribution at least 5 
of the Milky Way's satellites (the Large and Small Magellanic 



Clouds, Ursa Minor, Carina and Fornax) have proper motions con- 
sistent with "rotational support" of the so-called "disc-of-satellites" 
( |Metz et aT} 2008 1, an observation recently attributed to the angular 
momentum bias of satellites falling in along their hosts main axis 
l |Deaso n et al. 201 1). Although examining a different environment 
and mass range of haloes, [Tormen| < [1997[ > was the first to suggest 
that galaxies fall into clusters anisotropically. Further to their work, 
the idea that satellites are preferentially (not randomly) accreted 
and retain a memory of the large scales from which they came has 
also recently been suggested ( K nebe et al.||2004| |Libeskind et al.| 
[20TT][\faa^Ciro et al.|20l"T> . 

On larger scales a number of observational studies have hinted 
that galaxies separated by as much as ~ 30 Mpc (or 100s of times 
their virial radii) may also exhibit alignments albeit with the large 
scale distribution of matter (Hirata & Seljak 2004, Hirata et al. 
[20071 pandelbaum et aT1|2006| |Okumura et al.||2009> . On these 
scales, matter clusters into a "cosmic web" defined either geomet- 
rically (e.g. Sousbie et al.|2008) l or dynamically (e.g. Ha hn et al.| 
[20071 |Forero-Romero et al 12009) . 

Yet what lacks from the above studies is a clear link between 
the large scale cosmic web, the orientation of central galaxies and 
the preferential bias in satellite entry points, angular momenta and 
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Figure 1. A 0.25/i -1 Mpc slice through the simulations box of side length 
6Ah~ 1 Mpc. In the left panel, we show the mesh cells associated with knots 
(black), filaments (light gray) and sheets (dark gray). On the right we show 
the distribution of DM haloes in the same simulation slice. 

z = position. Some studies, such as Ara gon-Calvo et aL]^2007} 
and Hahn et al. ( 2007 ) have used methods based on the tidal field 
or the Hessian of the density field to correlate intrinsic halo proper- 
ties with the cosmic web. | Aragon-Calvo et a l. ( 2007 ) found a mass 
dependence for the (perpendicular) correlation between halo spins 
and cosmic structure. Hahn et al. (2007 1 examined how any correla- 
tion varied with redshift, extending their work in |Hahn et al.| ( |2010"| >, 
where they attempted to establish a link between central galactic 
discs and the cosmic web. In that paper they found an alignment 
between the spin of gaseous discs and the intermediate axis of the 
of the large scale tidal field, consistent with linear tidal torque the- 
ory (e.g. |Navarro et al.|2004|[Lee & Erdogdu 2007 1. However none 
of these studies examined the orbits of satellite galaxies. 

In this Letter we extend and complement the work of Hahn et 
al. (2007, 2010) and |Aragon-Calvo et al."|l2007| > in two ways. First 
we examining the orientation of the halo with respect to the velocity 
shear tensor instead of the tidal tensor. This method allows us to 
obtain much higher spatial resolution for our web classification. 
Secondly, we look at the orientation with respect to the cosmic web 
of substructure angular momentum, in a bid to understand how the 
internal dynamics of haloes and the substructures resident within 
them are effected by and reflective of the large scale cosmic web. 



2 METHODS 

We use a DM only Af-body cosmological simulation run assum- 
ing the standard ACDM concordance cosmology (e.g. WMAP5, 
|Komatsu et al.| |2009). These assume a flat universe with cos- 
mological constant density parameter J1a = 0.72, matter den- 
sity parameter f2 m = 0.28, a Hubble constant parameterized by 
Ho = WOh km s _1 Mpc^ 1 (with h = 0.7), a spectral index of 
primordial density fluctuations given by n s — 0.96, and mass fluc- 
tuations given by as — 0.817. 

The simulations span a box of side length 64/i _1 Mpc with 
1024 3 particles, achieving a mass resolution of ~ 1.89 x 
10 7 7i _1 Mq and a spatial resolution of l/i _1 kpc. The publicly 
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Figure 2. Normalized histograms of the cosine of the angle between the 
spin axis of the halo (Ldm. left), the orbital angular momentum of sub- 
structures ( J su b right) and the principal axes of the cosmic web in sheets 
(top), filaments (middle) and in knots (bottom). We show the angle formed 
with the long (ei, black), intermediate (e>2, red) and short (S3, blue) axes of 
the cosmic web. A uniform distribution of angles is shown by the flat line 
atP(|cos0|) = 1. 



available halo finder AHF (Knollman n~& Knebe||2009"l > is run on 
the particle distribution to obtain a halo catalogue. AHF identifies 
haloes and subhaloes in the simulation by searching the particle 
distribution for local density by maxima and checking that parti- 
cles within the virial radius are gravitationally bound to the host 
structure. Substructures are identified as haloes whose center is lo- 
cated within the virial radius of a more massive parent halo. 

A grid based velocity field is constructed by the "clouds in 
cells" (CIC) algorithm. The velocity of each cell is computed by 
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summing the momentum of all the particles in each cell and then 
dividing by the total mass in that cell. A normalized shear tensor 
is defined at each grid cell point by T, a = — jjj^ (§^T + 5^")' 
where a and /3 are the x, y, and z components. The Ho normaliza- 
tion is used to make the tensor dimensionless and the minus sign is 
introduced to make the positive eigenvectors correspond to a con- 
verging flow. The shear tensor is diagonalized at each grid point to 
obtain the eigenvalues Ai and the corresponding eigenvectors &i, 
(where i = 1, 2, 3 and Ai > A2 > A3). The web classification at 
each grid cell is done by counting how many eigenvalues are above 
a threshold A t h = 0.07, (chosen in Hoffman et al 201 1 to visually 
reproduce the cosmic web). If 0, 1, 2, or 3 eigenvalues are above 
Ath we classify the grid cell as a void, sheet, filament or knot re- 
spectively. The CIC velocity field is then gridded into 256 3 cells, 
each with a length of 0.25 h~ 1 Mpc and is then Gaussian smoothed 
with a kernel of 0.25 /i _1 Mpc so as to suppress numerical arti- 
facts and spurious preferred directions introduced by the cartesian 
grid. Note that the CIC algorithm inherently smoothes the discrete 
particle distribution across eight cells, two for each of the three 
dimensions. This intrinsic smoothing allows us to use a gaussian 
smoothing of 1 grid cell while still suppressing artificial grid ef- 
fects. 

Note that the number of positive eigenvalues corresponds to 
the number of eigenvectors along which mass is moving "inwards" 
or collapsing, while the number of negative eigenvalues indicates 
expansion. In this work we assume that if < Ai < A t h the eigen- 
vector is collapsing so slowly as to be considered expanding. Thus 
knots and voids are collapsing and expanding in all three directions 
respectively. Filaments, expand along their long axis, £3, while col- 
lapsing along ei and while sheets collapse along their short axis 
ei, and expand in the plane defined by 62 and £3. 

We focus on the orientation of halo angular momentum (Ldm) 
and the orbital angular momentum of subhaloes ( J su b) with respect 
to the cosmic web. Often the word "orientation" is used to mean 
the orientation of these vectors. For the study of Ldm only virial- 
ized haloes more massive that 1O 95 M are considered; this returns 
~ 10 6 haloes. Of these, ~ 2.7 x 10 4 haloes contain ~ 1.2 x 10 5 
substructures more massive than 30 particles, that are then used 
in the study of J su b- Only subhaloes with >30 particles are used 
here (AHF returns converged results for haloes above this limit, 
see |Knebe et al.|2011[ l, but note that increasing or decreasing our 
subhalo resolution to 20 or 200 particles, has no effect on our re- 
sults. 

The alignment of Ldm and J su b is measured with respect to 
the eigenvectors of the shear tensor. To each host halo a cosmic web 
classification can be assigned by finding within which grid cell the 
halo is located. It is the orientation of halo properties with respect 
to that grid cell's orfhonormal basis of the V-web that is examined 
here. 

For the purpose of this Letter we examine only haloes that 
are found at z — within sheets, filaments and knots. We ignore 
the orientation of haloes with respect to voids for two reasons (see 
table[T]and discussion in section[3]for details): the number density 
of haloes in voids is exceedingly small when compared to the other 
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Knots 


Filaments 


Sheets 


Voids 


Vff 


0.5% 


4.5% 


26% 


69% 


Nf 


6% 


28% 


45% 


21% 


M f 


25% 


44% 


26% 


5% 


N 


17 


6.1 


1.7 


0.3 


P 


72 


9.4 


0.96 


0.08 



Table 1. Properties of haloes and the cosmic web. We present: the volume 
filling fraction Vt t (the fraction of the total volume occupied by each web 
type), the number fraction Ne (the fraction of haloes assigned to a given 
web type), the halo mass fraction Mt (the fraction of all halo mass, in 
haloes of each web type), the number density of haloes TV (the number of 
haloes in a given web type divided by the volume occupied by that web type, 
normalized by the mean number density of haloes), and the mass density of 
haloes p (the total mass of haloes in each web type divided by the volume 
occupied by that web type, normalized by the mean halo mass density). 

web classifications and the haloes in voids tend to be low mass - 
little of the mass bound to haloes at z = is in voids. 



3 RESULTS 

Fig. [I] shows a 0.25 h^ 1 Mpc slice through the centre of our sim- 
ulation box. Since each point (grid cell) in space is given a web 
classification, we can visualize the cosmic web by dividing it into 
its constituent components. Fig. [T] reveals the nature of the cos- 
mic web: knots are fragmented over dense clumps in the puffier 
filaments, which themselves inhabit the centers of even less dense 
sheets. The filamentary appearance of sheets in Fig. [T] is a result 
of the planar cut of the two dimensional manifold. Although it ap- 
pears that sheets dominate the environment of haloes this is only 
due to their large volume. In Table [T] we present the properties of 
the cosmic web in terms of various filling factors. 

This table reveals the nature of the halo distribution - for ex- 
ample, although voids occupy 69% of the simulation volume, only 
21% of the haloes are found in voids. Furthermore, these haloes 
tend to be small, making up just 5% of the mass bound to haloes. 
Finally since the volume occupied by voids is so large, the number 
density of haloes here is 0.3 times the mean while the mass density 
is just 0.08 times the mean mass density of the haloes in the uni- 
verse. Conversely for knots, we see that although they occupy just 
0.5% of the simulation's volume, they contain a quarter of the mass 
bound to haloes and are 72 times denser then the mean density. 

Fig. [2] shows the normalized (differential) distribution of the 
angles formed between the principal axes of the cosmic web with 
the halo and satellite orbital angular momenta Ldm and J su b for 
haloes in a variety of different cosmic web environments. We per- 
formed Kolmogorov-Smirnoff (KS) tests to check the null hypoth- 
esis that the distributions are consistent with being drawn from a 
random (uniform) distribution. With the exception of the Knots, 
the null hypothesis can be ruled out at high (>99.99%) confidence 
levels. We examine these panels in detail below. 

In Fig. |2jt-c, we show the alignment of Ldm with the three 
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principle axes of the cosmic web. In sheets (Fig. [2^) we see a clear 
alignment with e$ and a (slightly weaker) perpendicular alignment 
with ei and e%. This indicates that a halo's spin axis lies in the sheet 
along the axis that is collapsing slowest (£3) and perpendicular to 
the sheet's short (normal) axis (ei). KS tests rule out uniformity 
at a high confidence level for all three axes. When we examine 
the orientation of Ldm in filaments (Fig. |2|)) we find a general 
weakening of the alignments found in sheets, but the same over 
all picture prevails - the halo spin axis points along the long axis 
of the filament (£3). Again, KS tests rule out uniformity at a high 
confidence level. The alignment weakens further in knots (Fig.|2j;) 
where the distribution of Ldm is fully consistent with a random 
orientation with respect to all three principal axes (KS probabilities 
of ~ 20%). 

Haloes that live in sheets are aligned perpendicular to the 
sheet's short axis, while filament haloes are aligned parallel to the 
filamentary axis, a result consistent with previous work (e.g. |Hahn| 
|et al.|20 07, Aragon-Ca lvo et al.|2007| >. In sheets, this is because the 
short axis is the axis of collapse, hence the axis from which mat- 
ter is being gravitationally pulled along. If the sheet is squeezed 
along this axis, it follows that the angular momentum will be per- 
pendicular to ei. Yet this interpretation weakens in filaments (as 
there are two "short axes") and disappears in knots. As haloes are 
drawn into filaments and then knots, their angular momentum be- 
comes randomized, most likely due to the highly non linear chaotic 
motions present in denser environments. 

Fig.|2jl-f shows the alignment of J su b, the orbital angular mo- 
mentum of substructures with the cosmic web. It is immediately 
clear that that J su b is always aligned parallel to the intermediate 
axis of the cosmic web - this alignment is strongest in filaments 
yet seen in all web classifications at significant KS confidence lev- 
els and is indicative of similar findings by Navarro et al. ( 2004 1. In 
both sheets and filaments however, there is additionally a (slightly 
weaker) perpendicular alignment with the small and long axes of 
the web. This indicates that substructures are accreted (and are thus 
moving) in a plane fairly well defined by ei (the sheet's short axis) 
and £3 (the sheet/filaments long axis). 

Remarkably the perpendicular alignment with £3 radically and 
significantly flips in knots (Fig. |2]). This indicates that although 
substructures are accreted in the ei-iz plane, by the time they land 
in knots the plane has nutated such that J sat is perpendicular to just 
e\, the main axis of the knot. The plane of accretion in filaments 
is well defined - the e,\-ez, while in knots this plane is allowed to 
nutate about e2- 

It is important to note that the significant perpendicular align- 
ment of J S ub with ei in knots is strong evidence that the cosmic 
web is well defined in knots. This indicates that the uniform distri- 
bution of angles seen between the cosmic web and Ldm in knots 
(Fig.[2j;) is indeed due to a randomization of Ldm. 

We close this section by noting that we also examined the ori- 
entation of substructure position with the cosmic web. Since the 
position vector is by construction perpendicular to J au b we find a 
mirror image panels (d)-(e), in otherwords, the positions of sub- 



haloes tend to align with <k\ and 63 in sheets and filaments and just 
with ii in knots. 



4 DISCUSSION AND CONCLUSIONS 

We have examined the orientation of DM spin axis (Ldm) and the 
orbital angular momentum of substructures (J su b) with respect to 
the cosmic web in sheets, filaments and voids. Our main findings 
arc 

• Substructures orbit such that their angular momenta tend to be 
aligned with the intermediate axis of the web structure they are in. 
If the second axis of the web is parallel to the intermediate axis 
of the material that collapsed to form the halo at turnaround, then 
our result is consistent with tidal torque theory (e.g. Nava rro et al.| 
|2004| >. This occurs because the velocity shear is closely related to 
the initial conditions. 

• In sheets and filaments J su b tends to be perpendicular to the 
short (ei) and long (£3) axis of the cosmic web, while in knots it is 
perpendicular only to the long axis. This indicates that in sheets and 
filaments, substructures move in just one plane (the ei-e3) while in 
knots the plane of their orbit is allowed any configuration parallel 
to the long axis, ei. This transition is most likely due to the ran- 
domization of orbits in the highly non linear halo interiors. 

• The spin axis of a DM halo tends to align with the short axis of 
the cosmic web in sheets and in filaments, but is consistent with a 
random orientation in haloes associated to knots. Since in sheets the 
spin axis lies in the plane, and assuming that in sheets the long and 
intermediate axes are degenerate, this result is also consistent with 
|Navarro et aL]pt304) ( given that there is some correlation between 
galaxy and halo spin). 

Our results indicate that the large scale velocity shear tensor is 
still a determining factor in the orientation of the angular momenta 
of substructures. Although Knebe et al. (2004) already claimed that 
the position of a subhalo at apocenter recalled its filamentary accre- 
tion even after several orbits, this link was not quantitatively estab- 
lished by them. This result is also consistent with a number of re- 
cent studies (Libeskind et al. 201 1 , Vera-Ciro et al. 201 1 1 that have 
argued that a memory of the large scale environment is imprinted 
on z = subhalo orbits and DM haloes. Yet none of these stud- 
ies were able to draw a direct parallel between the large scales and 
satellite orbits at z = 0, instead they relied on tying z — prop- 
erties of satellite galaxies to the merger history of the host halo. 
Indeed our result is surprising since we do no such tracking back in 
time. Instead our result indicates that the large scale velocity field 
is reflected in the orbits of subhaloes at z = 0. 

Since a subhaloes z — position is perpendicular to its an- 
gular momentum, our findings help understand the relationship be- 
tween the anisotropic z — spatial distribution of satellite galaxy 
in simulations (Knebe et al. 2004; Libeski nd et al.|200 5, Zentner 
et al.|2 005 , Libeskind et a H2007[> and observations (e.g. |Sales & 
Lambas 2004, Kroup a et al.||2005] a mong others) with their ap 



parent coherent motion ( |Metz et al.|[2 008 2009 ; Libeskin d~et al.| 
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[20091 IDeason et afpOn) . Subhaloes are assumed to host lumi- 
nous satellite galaxies that do not enter the halo randomly but in 
some correlated way (e.g. |Tormen| 1997j |Li & Helmi|2008l |Knebe| 
|et al.|20PT) . Our findings indicate that the orientation of the non- 
randomly accreted subhaloes is determined by the cosmic web at 
z = 0, thus connecting the large scale sheets, filaments and knots 
with the sub-virial radius scale of subhalo orbits. In summary, The 
non-random infall points of subhaloes is both reflective of the large 
scale structure and frozen into the z = angular momentum of a 
subhalo, while, due to long dynamical times, the 2 = positions 
of satellites are susceptible to randomization processes. 

The strength of the alignment is likely to depend on a number 
of factors - mass of halo, redshift, and how old a given parent halo 
is. A comprehensive study of the nature of haloes in the cosmic 
web is forthcoming (Metuki et al in preparation). 
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